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Abstract. We present a detailed study of the spherically symmetric solutions in 
Lorentz breaking massive gravity. There is an undetermined function 11 ) 1 , 11 ) 2 , 11 ) 3 ) 
in the action of Stiickelberg fields = A 4 / d^Xy/~^gT, which should be resolved 
through physical means. In the general relativity, the spherically symmetric solution 
to the Einstein equation is a benchmark and its massive deformation also play a 
crucial role in Lorentz breaking massive gravity. T will satisfy the constraint equation 
Tq 1 = 0 from the spherically symmetric Einstein tensor Gj = 0, if we maintain that 
any reasonable physical theory should possess the spherically symmetric solutions. 
The Stiickelberg field q i l is taken as a ’hedgehog’ configuration </>* = (j)(r)x l /r , whose 
stability is guaranteed by the topological one. Under this ansatz, Tq 1 = 0 is reduced 
to dT = 0. The functions T for dT = 0 form a commutative ring . We 
obtain an expression of solution to the functional differential equation with spherically 
symmetry if T £ . If T £ R^ and dT/dX = 0, the functions T form a 

subring S^ C R^. We show that the metric is Schwarzschild, Schwarzschild-AdS 
or Schwarzschild-dS if T £ SW When T € R^ but T ^ S■ F , we will obtain some 
new metric solutions, including the furry black hole and beyond. Using the general 
formula and the basic property of function ring R^, we give some analytical examples 
and their phenomenological applications. Furthermore, we also discuss the stability of 
gravitational field by the analysis of Komar integral and the results of QNMs. 


PACS numbers: 04.50 Kd, 04.20 -q 

1. Introduction 

It is an interesting question whether general relativity (GR) is a solitary theory from 
both theoretical and phenomenological sides. One of modifying gravity theories is the 
massive deformation of GR. A comprehensive review of massive gravity can be found in 
|T]. We can divide the massive gravity theories into two varieties: Lorentz invariant type 
(LI) and Lorentz breaking type (LB). LB massive gravity can be physically viable, which 
is free from pathologies such as the ghosts, the van Dam-Veltman-Zakharov (vDVZ) 
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discontinuity due to extra graviton polarization and strong coupling at the low energy 
scale (see e. g., Refs. EE] and references therein). In most of LB cases the energy 
cutoff ( mM p i )5 is quite high in comparison with that in LI cases, such that LB massive 
gravity can avoid some phenomenological shortcomings from LI massive gravity [4]. 
For instance, there is a class of LB massive gravity models in which the solar system 
constraints are satisfied for larger graviton mass [5]. The relic of gravitational waves 
engendered during inflation epoch may constitute today the cold dark matter in the 
Universe and would bring to a distinguishingly monochromatic signal in the gravitational 
wave detectors [6]. However, the significant contribution of massive gravitons to the local 
dark halo density have ruled out essentially [7]. Recently, the massive gravity theories 
with five propagating degrees of freedom were made in a series of paper [SJ \9\ [10j • The 
Lorentz breaking in their gravitational sector is inescapable to concede a nontrivial 
cosmology for a spatially flat Universe [3]. 

In GR, the spherically symmetric vacuum solution to the Einstein equation is a 
benchmark and its massive deformation also play a crucial role in LB massive gravity. 
The exact spherically symmetric solutions in some massive gravity models show a variety 
of different features mmmm- They differ from the Schwarzschild solution, 
which depend on two parameters: the mass M and the scalar charge S. Especially, 
the solutions show a nonanalytic hair in the form of a powerlike term r - A m- The 
solutions may possess a horizon depending upon the parameters M and S, which become 
candidates of black holes. 

In most theories of LB massive gravity, the Lorentz group is broken down to a SO (3) 
rotation group. A class of promising LB massive gravity theories is especially interesting 
where the ghosts are absent because there is the residual symmetry x 1 —> x l + f l (t) on 
flat space for three arbitrary functions ff{t). In the picture of Stiickelberg fields 0°, 0* 
(i = 1,2,3), this corresponds with the internal symmetry 0* —y 0* + £*(0). Thus, the 
mass term should be a function T of the dimensionless quantities X = g^ u d ^0°c00° 
and W = ( g^ — X~ l d^<f 0 d v <fP)d^(() l d u <fP [2]. The authors of Refs. [Ill ITT] studied the 
black hole solution of models described by the designated function T. 

In this work, we release from the limitations of designated function T so that the 
Einstein equations become the functional differential equations at present. We have 
found the functional formula of spherically symmetric solutions. In GR, the Birkhoff’s 
theorem is the statement that the Schwarzschild metric is the unique vacuum solution 
with spherical symmetry. The standard Birkhoff’s theorem cannot be applied because 
there are the Stiickelberg fields in the massive gravity. Therefore, the function T must be 
not arbitrary, but restrictive if we consider a static spherically symmetric system. The 
Stiickelberg held 0* is taken as a ’hedgehog’ configuration, whose topological property 
guarantees that the configuration is stable. The proper functions T should satisfy 
dT = 0 from the condition of spherically symmetric solution Tf = 0 where Tf is the 
energy-momentum tensor of the Stiickelberg fields. Under the hedgehog ansatz, the 
functions J- form a commutative ring R?. The basic property of R T affords us how to 
choose T for the various solutions. Before leading to new solutions, we first checked a 
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well-known example using our universal formula of the spherically symmetric solutions. 
The check has a double significance, which not only checked the known solution but also 
checked the formula itself. If the functions T G R T and dTjdX = 0, T form a subring 
S T C R T . We find that the metric is Schwarzschild, Sch-AdS or Sch-dS if X G S T . We 
obtain new analytical solutions as well as the well-known example when T G R T but 
X . Especially, one of new solutions has interesting astrophysics applications. 

The paper is organized as follows: in section 1 we introduce the Lorentz breaking 
massive gravity. In section 2 we shall answer how to choose X and show that the proper 
functions X form a commutative ring Vft. Furthermore, we shall give the functional 
differential equations and the expression of static spherically symmetric solutions. Then, 
we shall give some analytical solutions, classification of the functions X, and discuss 
stability of the Stuckelberg held configuration in section 3, analyse the stability of 
gravitational held of these solutions in section 4 and study the phenomenological 
consequences for new solution in section 5. Finally, section 6 is devoted to the conclusion 
and discussion. 

2. Static spherically symmetric equations and solutions 

In this work, we discussed the spherically symmetric solutions in a wide class of 
promising LB massive gravity theories, where the only propagating degree of freedom is 
a single graviton [[3j. Following Stiickelberg [15], general covariance is restored by a set 
of Goldstone helds, which are nonlinearly interacting with the metric held. The action 
can be written as S = Sgr + S^, and 

= A 4 J ftxft^x(x, w ij ). (l) 

describes a set of four Stuckelberg helds ft, ft (i = 1,2,3) transforming under 
a cliffeomorphism (cliff) 8x M = ft(x) as simple scalars. One can structure cliff 
invariant function X using dimensionless quantities X = g^ J,u d fl ftd l/ ft and lf IJ = 
(g^ u — X~ 1 d^ftd u ft)d fl ftd u ft. The helds ft and ft have spacetime dependent vacuum 
expectation values 

g,w = = at i ft = bx *■ ( 2 ) 

If a function T is chosen, we can redefine helds ft and ft so that a = b = 1. The 
action (JT]) is invariant under ft —> ft + ft (ft), where ft are arbitrary functions of 
ft. This symmetry turns out that the perturbations about the vacuum solution (J2]) 
are nonpathological [3j. One requires that the background breaks only the invariance 
under boosts but not under the whole Lorentz group. Then, the function T depends 
on W through three combinations w n = TrW n (n = 1, 2, 3) which is necessary for the 
background to be invariant under the SO (ft) symmetry in the ft internal space [Tl[ 14] , 
The constant A has the dimension of mass, and S ( p is understood as the low-energy 
effective action valid below the scale A. 
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2.1. How do we choose X? 


The static spherically symmetric ansatz can be written as follows m 
ds 2 = aft)dt 2 — ftft)dr 2 — r 2 dfl 2 , 

ft — t + hft), 

ft = -ftftxft (3) 

r 

where ft is the so-called ’hedgehog’ configuration. We note that the ansatz (J2D is not 
the most general one consistent with the symmetries. The spatial Stuckelberg field is 
ft = bx l in the end (see Lemma 1 in Appendix A). In other words, they are actually 
fixed to their unitary values. However it may be that the static spherically symmetric 
ansatz is written in different coordinates in unitary gauge [16]. 

This configuration contains two additional radial functions, lift) and cf>ft ) as 
compared with GR. Since lift) = constant corresponds to vacuum solution, we might 
as well assume h'ft) ft 0 for nontrivial solutions. Using ([3]), we can express w n and X 
by the following relations 

W\ — — ft 2 {aft)~ l X~ 4 — 2 ftr~ 2 , 
w 2 = ft 4 (aft)~ 2 X~ 2 + 2 ftr~ 4 , 

w 3 = - ftftaft)~ 3 X~ 3 - 20 6 r" 6 , (4) 


and 


X = 


ft — ah' 2 
aft 


(5) 


The non-zero components of energy-momentum tensor Tf for the Goldstone fields are 
given by the following expressions 

T 0 ° = A 4 [X + 2(aftXft 1 ft 2 X 1 - A(aftX)~ 2 ft A X 2 + 6(aftX)~ 3 ft e X 3 \, (6) 

2 A 4 A' 

ft = - — [X x + {aft)- 1 X~ 2 ft 2 X l - 2 (aft)- 2 X~ 3 ft 4 X 2 + 3(aft)- 3 X~ 4 ft e X 3 ], 


ft 


ft = \ 4 [X - 2X x X + a~\2X x + 2{aft)~ l X~ 2 ft 2 X - 4 (aft)- 2 X~ 3 ft 4 X 2 
+ 6(aft)- 3 X~ 4 ft 6 X 3 )\, 

9 o ft ft dft 

Tl = T 3 = X + 2ftX - 4ftX 2 + 6 ^X 3 , 

where Xi = dX/dwi and Xx = dX/dX. The non-zero components of the Einstein 
tensor are 

I r 

,2 


(7) 

( 8 ) 
(9) 




r x_ 1 n _ a + raft 
T1 r 2 aft J 


( 10 ) 

( 11 ) 


1 r a'+ ra" .2 a + raft.. 


_ 

(jr 2 ^3 


4r aft 


( 12 ) 
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From the Einstein equations Gf = M p 2 Tf and Gj = 0 we have Tf = 0. It is easy to 
prove that T 0 ° — Tf is proportional to Tf and T 0 ° = Tf implies Gg = G\. Thus, we have 
a(r)/3(r ) = 1, and the Einstein equations can be rewritten as 


V 

-h 

r 
, // 


, / 


(X-2XXx) = 0, 


- , - 4 0 4 T 60 6 
Y + 7 - Y ( ^ + ~ + - 0) 


i/2 


i/4 


i/6 


^ + X ^ 1 _ 2 ^ P ' F ‘ 2 + ^ Y 1 ^ 3 - °' 

From aft = 1, (j4j) and ([5]) can be reduced to 

wi = -( 0 , 2 X -1 + 2 0 2 r~ 2 ), 

w 2 — (0 /2 X -1 ) 2 + 2(0 2 r -2 ) 2 , 
m 3 = -[(0 /2 X” 1 ) 3 + 2(0 2 r- 2 ) 3 ], 


(13) 

(14) 

(15) 


and 


X = 


a 


l h ' 2 


a 


(16) 


(17) 


Eqs. (fT3]h(fT5l) are the system containing three differential equations if and only if 
we take a fixed function T. From this point of view, the system contains three equations 
for three unknowns, a, <f> and X. Even though it is a highly non-linear coupling system, 
it is solvable. Our purpose is to find an expression of spherically symmetric solution and 
analytical solutions for some particular choices of the function T. We adopt a strategy 
to achieve this purpose as follows: The first step is that (fl5|l is reduced to an algebraic 
equation if we assume T is a polynomial of X -1 , w\, w 2 and w%. Thus, we show that 
there is only solution 0(r) = hr and (Il3l) - (fl5l) become a system containing two functional 
differential equations. Secondly, we find an expression of a(r) for any polynomial X. 
Thirdly, we obtain X using the relation between a and A". Finally, we obtain 0° using 
the relation between X, a and 0°. 

To solve the functional-differential equation (fT5lh we realize that only one of 
variables X,wi,w 2 and u> 3 is effective. This occurs if there are subsidiary relationships 
constraining all the variables to have values dependent on the value of A". These 
relationships may be represented by equations 


Wi = Wi(X),(i = 1,2,3) (18) 

There is only solution for (Il8l) . 

0(r) = br, (19) 

where b is an undetermined constant and it will be used in (fl6l) . Fortunately, we can 
prove that (1T91) is the only solution to Einstein equations (fT3lh(fT5l) under the hedgehog 
ansatz ([3]) and T is a polynomial of X -1 , wi, w 2 and w 3 (see Lemma 1 in Appendix A). 
Thus, (TT5h can be rewritten as 


dX = TxdX + Xidwi + X 2 dw 2 + X^dwz = 0 . 


( 20 ) 
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We note that the gravitational fields obey the Bianchi identity = 0 at present 

which implies the conservational requirement V^X 7 ^ = 0. The energy- moment urn tensor 
is really conservation in the 0 = br, so there is consistency between the Bianchi identity 
and the additional restriction (fT5l) on R. 

We consider R is a polynomial below, 

jr = a iohi2h(^y ow i w W3 ( 2i ) 

where Iq, ■ ■ ■, l 3 are positive integers. Mathematically, a commutative ring is a set 
with two laws of composition + and x, called addition and multiplication, that satisfies 
rules: (i) With the law of composition +, is an abelian group that we denote by Rf F+ \ 
its identity is denoted by 0; (ii) Multiplication is commutative and associative, and has 
an identity denoted by 1; (iii) With the law of distribution (P + f j )f k = f\f k + f J .f k and 
/ fc (/* + /‘ ? ) — f k f l +f k P f° r all f\ f J ! f k in R T . In our case, all solutions (l2lf of dR = 0 
form a function ring R T . Remarkably, there exist different solutions for a{r ) and h(r), 
even though 0*(r) is provided with same form. The algebraic feature of R^ shall aid us to 
find new solutions a(r) and h{r). Especially, all solutions T = ... ^ aoip 2 i 3 w l f w l f w l f 

form a subring S T C R T . If T G S T , we will prove that the solutions of Einstein 
equation are the Schwarzschild, Sch-AdS or Sch-dS respectively. 


2.2. Functional differential equations 

Under the hedgehog ansatz ([3]), there is the only solution 0 = br for Einstein equations 
(fT3lb (TT5l) . Thus we choose T G R T , the Einstein equations are reduced to the following 
functional differential equations 


= 0, 


Ji' 

— + 

Y + 7-T^> = °' 


( 22 ) 

(23) 


and a/3 = 1. Here m 2 = 2A 4 /M 2 l and the functionals K/(/F) and R{T) are dehned as 
JC(R) = R - 2XR x , (24) 

H(R) = R + 2 b 2 R 1 - 4 PR 2 + Qb 6 R 3 . (25) 

f)22|) and (123|) are well dehned. Since R is a function of X, w 1 , W 2 and uq, the functionals 
)C(R) and R(R) are also ones. Furthermore, w 1 , w 2 and w 3 are the polynomials of X~ x 
from (fTBjl and (flUD . so /C and FL are also looked upon as functions of X. Thus, A" can 
be expressed as a function of a, a' and r from (|22|) . and (l23|) becomes an ordinary 
differential equation of a. Therefore, we can obtain a from (125]) . then we have the 
solution 0° from (fT7|) and (122|) . 

We solve formally (]22l) as follows 

A = K-\y), 


(26) 
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where 


2 ,a' 

?/ = -?(- + 


a 


m 2 ' r 


(27) 


According to the fundamental theorem of algebra, every nonconstant polynomial has a 
complex root. However, if T £ R T , T + c is also in R T so that we always obtain real 
Kr x {y). This result can be generalized for any differentiable function in R T . 

There is a mathematical identity relation 


2a" + 4r W = m 2 ry' + 2m 2 y , 


(28) 


which is the key to the solving process. Using (1281) . (123D can be rewritten as 

r -y' + y-n(JC~ 1 (y)) = 0, (29) 

where TL{1C~ x {jj)) is also the functional of T, and we can regard it as a function of y. 
Therefore, ()29li is a separable first-order equation of y in reality. When we solve (|29jl . 
there are two cases: (i) TL = /C; (ii)'H ^ 1C. We have Tx = 0 if 'H(-A) = tC(J-) (see 
Lemma 2 in Appendix A). 


2.3. The Schwarzschild, AdS and dS solutions 


If T G S T , we have TL = 1C. In this case, (129]( has only constant solutions, y = const., 
which correspond the Schwarzschild, Sch-AdS or Sch-dS solutions respectively from (127)1 . 
That is to say 


a = 



for y = 0, 
for y = 


From ([5]) for the Schwarzschild solution, we have 

, , , „ (2 -/C _1 (0))r s , /(1 ^ /C- 1 (0 ):r s ^ 

A + r s cosh 1 B + -—-In (A + Cr H-——) + u 0 


2 C 


A + r 9 \n CL , V ^ ) + Uo 


2 C 


h(r) = r + r s In (r - r s ) + u 0 


'r + y/r s 

where w 0 is an integral constant and 


A — [(l-K. 1 (0))r 2 + /C 1 (0)r s r , ]3, 
(2-/C- 1 (0))r s r + /C^ 1 (0)r 2 
|/C-i(0)|r s (r-r s ) 

C = [l-/C- 1 (0)]5. 


(30) 

for /C-^O) < 1 
for /C -1 (0) = 0 (3i) 
for /C- 1 (0) = 1 


It is worth to note that 0°(r) is not real solution if LC x (0) > 1. 


( 32 ) 
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In the case of dS solution a = 1 — fi 2 r 2 , we require 1C x (—^-) < 1 for 0°(r) being 
real solution. From (J5]), we have 


( _ ^ ln ^ 1 “ h 2 ' r2 ) + u 0 


u o 


h{r ) = < 


-1 In [ A ^~ B + C \ C - D W^A-^h + 

2/i [ B(l + A + CyC + D J 1 

1 , A 
2p 1U B +U ° 

1 [ ln T B a + S - 2 \l Sin E] + u 0 for K~\-K) < 0 


for = 1 

v / 

for 0 < IC-y-K) < 1 
for = 0 


(33) 


2/i 5(1 + A + C) v m “ 

where «o is an integral constant and 

A = 1 + fir, 

5 = 1 — fir, 

C=[l-C- 1 (-^)(l-^r 2 )]i, 

rrr 


" m z 


E = 


k-h-K) 

\ rri, z / 


-/ir. 


(34) 


For the case of AdS solution a = 1 + /i 2 r 2 , 0°(r) is not real as /C _1 (^ 2 -) > 0. From 


), we obtain 


h(r) = < 


( 1 _! 

— tan /rr + Uo 

fi 


for Kr'iK) =0 


1 nr 

— tan 1 ; + 

t-i yi - A 


-ic-'iH 




1 n(\llC-^^f)(A - 1 )-B)+ u 0 for K,~\^) < 0 


(35) 


m- 


where w 0 is an integral constant and 
A = 


(36) 


Similarly, we can obtain h(r) for Sch-AdS or Sch-dS from (151) . In these cases, the 
metric is independent with the integral constant Uq, so it only represents a residual 
freedom for the Stiickelberg held 0°. 


2 . 4 . The expression of solution in the case of Tty 1C 

If T G R T but T 0 S T , we have Ti ^ 1C. In this case, we obtain the formal integration 
from (129]) as follows 

dy 

^(^ -1 ( 2 /)) -y 



( 37 ) 
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Thus we expressed y as a function of r using the corresponding relation (1HTD between r 
and y , i.e., y = y(r,u 0 ) where u 0 is an integral constant. (T27|) can be rewritten as 


a' a — 1 m 2 
r + r 2 = 2 

(38) 

then we have the formula of ordinary solution 


77^2 /* 

a = lH-/ r 2 y(r,Uo)dr. 

2 r ./ 

(39) 

From (j3J), (122|) and a/3 = 1, we obtain 


d(r) = ± [ —[1 — a/C _1 (//)]2 . 

.7 a 

(40) 


Thus, we finally obtain the expression of solution for the functional differential equations 
in the TL ^ /C case, 

,2 r 


a = (3 1 = 1 + ^- / r 2 y(r,u 0 )dr , 


(jr — t ± 
f ) 1 = for*. 


dr 


a 


[1 - a/C 1 (y)] 2 , 


(41) 


3. Analytical solutions beyond known solutions in GR 

3.1. The solutions of furry black holes 

To find new solution beyond Schwarzschild, AdS and dS, we must take T G R T but 
T ^ S T . In the simplest case, we choose that T satisfies /C [T) = cH(IF), c G 1Z. Thus, 
we have 


?*(£ 1 ( 2 /)) = c y, 

and (1291) is reduced to 

V + (1 - c)y = 0. 


The ordinary solution of 


y = w 0 r 2(c d, 


is 


(42) 

(43) 

(44) 


where Uq is integral constant. Substituting (1421) and ilTTIi into (TUT) , we obtain the solution 


a = r i = i_ri_4, 

dr, 


= t± I —[1 - a/C ( y )] 2 , 

J a 

<\? = fex*, 

2 

where A = —2c, S' = — (^7 and r s is another integral constant. 


( 45 ) 
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Next, we reconsider a well-known example HQ to check up the formula (I45p . Taking 
the function 


h 2 

T = c 0 ( — + w i) + ci(wf - 3wiw 2 - 6b 4 w 1 + 2 w 3 - 12fe 6 ), 

YY 


(46) 


where Co and C\ are dimensionless constants. Two additional conditions should be 
enforced on these two constants: (i) Non-tachyonic condition Co — 6ci > 0; (ii) No ghost 
condition Co > 0. Thus, we have 

6cifr 4 


£ (y) = 


2 cnb 2 


and c = 


(47) 


y + 2 cob 2 c 0 

Substituting (1441) and (1471) into (145]) . we obtain again the solution in m by the markedly 
dissimilar way, 


a = f3 4 = 1 -- — 


S 




= bx\ (48) 

We note that the dimension of Uq is [Length] A+2 in the furry black hole (145j) . The 
parameter u 0 satisfies u 0 > —if we take c\ = —■§?[■ For this black hole, we have 


u 0 = -—^rr- if it has solar mass and horizon 


By the straightforward calculation, we find a more simple solution J-(X, W\, w 2 ) 
which satisfies = cR{fF) as follows 


T = c 0 (— + tci) + ci(w\ - w 2 + Ab 2 wi + 6b 4 ), 

YY 


and 


c = 


2ci b 2 

c 0 


(49) 


(50) 


where dimensionless constants c 0 and ci should satisfy the constraints c 0 + 4ci > 0 and 
c 0 > 0. 

Thus, we have 

2 Cob 2 


K (&) = ,n 72 ’ 

y + 2 cob 1 

and we finally obtain the solution 

a = r 1 = 1 - - - 4. 


(51) 


= t± 
7 = bx\ 


dr 

a 




(52) 


Using the same procedure, we take 
6 6 


T = c 0 ( 22— + w 3 ) + ci(wi + 3b 2 wi + 3 b 4 )(wl - w 2 + Ab 2 Wi + 66 4 ), (53) 
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and the analytical solution is 


a = f3 1 = 1 -- — 


S 


= t ±j *[l-a(l + 

= bx\ 


-Jf°_ 7 .-A- 2 Wi§ 

6 c n 6 6 J J ’ 


(54) 


where S = A = and u 0 , r s are integral constants. 


_ 4 c\b 2 

3 c 0 


There are certainly the infinite solutions with T G 7 Z T and /C( J-) = cH(J-). 
Furthermore, each fixed c can generate a subgroup Sff which satisfies fC(J-) = cK[T) 
if T G Sf (see Lemma 3 in Appendix A). For any T that satisfies = cH(T ), we 

have the same form for metric 

a = f3~ l = 1--- j. (55) 

This metric differ from the Schwarzschild solution by an additional powerlike term r -A . 
Such solution may possess an event horizon depending upon the parameters of the 
function J r , which become candidates of modified black hole. In other words, such black 
holes can be described by two physical parameters: Schwarzschild radius r s and scalar 
charges S, so that they are dubbed furry black holes 


3.2. An analytical solution beyond furry black holes 

If we want find the solutions beyond furry black holes, we have to choose J- G TZ T 
and /C(J 7 ) 7 ^ cH{J-). Using the expression of solution (j4T]h we are able to obtain new 
solution. For instance, we take J- as follows 


b 2 

T = c 0 (— + wi + 2b 2 ) + ci(wi + b> 2 ){wl -w 2 + 46 2 u>i + 6b 4 ), 

Ji 

where c 0 and Ci are dimensionless constants, and b 2 = (c 0 /2ci)^. Substituting 
and (J25]), we have 

1 


and 


K{?) = 

U{F) = 


a?X 

1 

a 2 X 2 


where the parameter a 2 = (ci/2cq) 2 . Thus, we obtain the expression /C l (y) by 

KT\y) = 4_. 
ay 

Similarly, we also obtain the expression 'H{fC~ l {y)) by (l58li and (l59li 
H{K-\y)) = a 2 y 2 . 

Thereupon, (l23l) can be reduced to 


(56) 
into 

(57) 

(58) 

a. 

(59) 

(60) 


-/ + H - a 2 !/ 2 = 0, 


(61) 
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which is provided with the solution 

1 

y = 


a 2 — vfpr 2 ' 


(62) 


Here, u$ is an integral constant. We have the solution for this system as follows 
?-i _ 1 -2 


a = 0 _i = 1 — e 2 - 25 — 

r i/2 mr 


3 In ^° r ~~ ^ 


Uor + a 


iO 


= t± [ -[l-a(l-^r 2 )]i, 
J a a A 

/ Co .1 i 
= ( 2 ^ *’ 


(63) 

where e 2 = ^4 and r s is another integral constant. (163|) gets a new type of spherically 
symmetric metric which has interesting properties. Let us discussed its behavior and 
astrophysical consequences in Section 5. 


3.3. Classification of the functions J 7 

The functions T G 'Rs F may be classified under three types according to metric form, 
that is to say, they are classified by the relation between /C(IF) and 'H(iF). In the special 
case, the solution of furry black hole is reduced to Schwarzschild solution (5 — 0), Sch- 
AdS and Sch-dS solutions (A = —2, S ^ 0). Three types are 

Type (i) : In the case of /C( J r ) = 'H(iF), there are only Schwarzschild, Sch-AdS 
and Sch-dS solutions. 

Type (ii) : In the case of IC(iF) = cK{fF\ there are the solutions with metric of 
furry black hole. 

Type (iii): In the case of ^ c'H(iF), there are the solutions beyond the metric 

of furry black hole. 


3-4- Topological stability of Stilckelberg field configuration 

There is no way to smoothly deform the hedgehog configuration into configuration 
(p l = 0 everywhere. Thus, it is topological stable. The simplest example is the ’t 
Hooft-Polyakov monopole which appears when SO(3) is broken to U( 1). The three- 
component Higgs field 0 * for the monopolc takes hedgehog configuration mm - The 
generalization of’t Hooft-Polyakov monopole is 50(10) GUT solution which appears 
when 50(10) is broken to 50(3) x 0(1) [19]. Mathematically, topological classes of 
solutions can be identified with classes of homotopic two-sphere 5 2 . This is expressed 
by the statement that n 2 [50(3)/0(l)] = n 1 [0(l)]/n 1 [50(3)] = Z/Z 2 for the ’t Hooft- 
Polyakov monopole. The tensor monopole has also been discussed in a LB theory BUJ- 
We have chosen the function T depending on IL !J through three combination 
W n = TrlP™ (n = 1,2,3) which will be necessary for the background to be invariant 
under the 50(3) symmetry in the 0* internal space [11., [Hj. In the case of 0* = 0, 
LB massive gravity will degenerate into GR, and Lorentz symmetry is recuperated. We 
note that the homomorphic relation is well known between 5L(2, C ) and the Lorentz 
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group. Therefore, the universal covering group is 577(2) x SU{2) for the Lorentz group. 
Since 


n 


r S77(2) x SU{ 2). 


= T\ l [SO{3)} = Z 2 , 


(64) 


S79(3) 

the Stiickelberg field conhguration ([3]) corresponding to a monopole exists, and there is 
no way to smoothly deform the monopole into a conhguration with Lorentz symmetry. 

We note that the Minkowski solution ([2]) is consistent with the monopole 
conhguration </>* = bx l . Therefore, the topological stability guarantees the existence 
of spherically symmetric solutions including the trivial one. We have to analyse the 
stability of a non-trivial spherically symmetric solution in the perturbed state in the 
next section. 


4. Stability of furry black holes 

Isolated furry black holes are nonsingular outside the event horizon and fully 
characterized by the parameters r s and S'. However, a real furry black hole must not 
be fully described by these basic parameters and it is invariably in the perturbed state. 
We would know something about the interaction of black hole with its astrophysical 
environment so that we have to analyse the stability of furry black holes in the perturbed 
state. 

There are three stages when a black hole is perturbed. First stage is the rapid 
response at very early time, on which the initial conditions have a great effect. Second 
stage is quasinormal ringing phase, whose characteristic oscillation frequencies and 
damping times depend strongly on the quasinormal modes (QNMs). The QNMs are 
determined completely by the parameters of system, so they would carry significant 
information about the background curvature of the intervening spacetime. Finally, there 
is a tail stage, which decays approximately as a power in time owing to backscattering 
off the spacetime curvature [21] , QNMs have been extensively studied for various black 
holes by different methods m including the string black holes [23], the brane-world 
models [54] and the acoustic black holes [25] - In order to comprehend which of these 
furry solutions are stable, and thus could exist in nature, we should consider their QNMs. 
Stability is assured when all QNMs are damped. 

4.1. Classification of furry black holes 

Note that the spacetime with metric (1551) is asymptotically flat which tends to vacuum 
at spatial infinity if only A > 0. Furthermore, we only consider the solutions that there 
exists the event horizon. If so, we can divide the metrics into following types, 

Type IA : A > 0, r s > 0 and S' > 0; (65) 

Type IB: A > 1, r s > 0, S < 0 and r s >- --(S' — AS 1 )*; (66) 

A — 1 

Type IC : 0 < A < 1, r s < 0, S' > 0 and r s >- --(S' — AS)*;(67) 

A — 1 
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Type IIA : A > l,r s < 0 and S > 0; ( 68 ) 

Type IIB : 0 < A < l,ry > 0 and S < 0. (69) 

By the analysis of Komar integral and the numerical results of QNMs, we argue that 
Type I is stable and Type II is unstable. Complcmentarily, there is an event horizon iff 
r s > r crit for the Type IB and IC, and the critical value r crit = — AS')*. 

4.2. Komar integral for furry black holes 

The total energy (or mass) is a tricky notion in the massive gravity or GR, because 

energy-momentum is a tensor than a vector, and energy-momentum tensor T /lL/ only 

describes the properties of matter, not those of gravitational field. Fortunately, one can 
still define a conserved total energy by the Komar integral [26] if stationary spacetime 
is equipped with a timelike Killing vector field The Komar integral associated with 
can be written as 

E = [ d 2 x „V / '£", (70) 

4?tG Jqy, 

where T is the spacelike hypersurface, ttF is the unit normal vector associated with £, 
the boundary dT is equipped with metric 7 and cr M is outward-pointing normal 
vector associated with dT. When the integral (1701) is convergent we can use it as the 
definition of total mass in all stationary asymptotically flat spacetime. 

For the furry metric (l55]h the boundary dT can be taken two-sphere typically, so 

that 

r )\^ > dx l dx^ = r 2 (d6 2 + sin 2 9d<f> 2 ) (71) 

and the normal vectors, normalized to = 1 and 070 ^ = — 1 , have nonzero 

components no = « 2 , = —a ~2 with other components vanishing. Thus, we have 

W» cr *'V /i £ 1 ' = — V 0 ^ 1 . The Killing vector is = (1, 0, 0, 0), so we have V 0 ^ 1 = g 00 Tl 0 l C° = 
Putting them all together, the Komar integral of a furry black hole is 
If 1 \o 

E = J dOdcfr 2 sin 0 a' = — (r 8 + p^y). (72) 

We note that the Komar integral (1721) still depends on the radial coordinate r, which 
demonstrates a remarkable characteristic of the furry black hole. Obviously, the Komar 
integral is only convergent if A > 1. If so, we can use (172]) as the definition of total 
energy (mass) when r tends to infinite. On the contrary, we can’t define total energy 
(mass) if 0 < A < 1, even though the Komar energy defined well at the spacial finite. 

Crossing f = (—the Komar integral would change sign so that the attractive 
gravitational field converts into replusion or vice versa. For Type I case, the Komar 
integral is always positive definite on the outside of horizon 77 . For Type II case, the 
Komar integral is positive definite in 77 < r < f and negative definite in r > r. 
It is known rigorously that the Schwarzschild black hole is stable against external 
perturbation [57]. The results for QNMs support this conclusion, as shown by the 
negative imaginary parts of QNMs [28].The Schwarzschild black hole is a special case of 
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Type IA whose Komar integral is positive definite. The Komar integral is also positive 
definite outside horizon for Type IB and IC. However, the Komar integral is negative 
definite for Type IIA and Type 11B in the region r > r, which is an essential distinction 
between Type 1 and Type II. Thus, we speculate that Type I is stable and Type II is 
unstable which can be verified by the results of QNMs for various furry black holes. 


4.3. Quasinormal modes of furry black holes 


The QNMs for spin 0 and 1 may be derived assuming the coupling between the spin 0 
and 1 fields to gravity is the same as those in GR. Since Lorentz invariance has been 
broken by mass term and therefore the calculating QNMs is not obvious in the case of 
gravitational perturbations. The perturbation equation for a furry black hole can be 
reduced to Schrodinger-like form for stationary background 

AAQ + [„* - V(rMr.) = 0, (73) 

where the tortoise coordinate r* maps the semi-infinite region from the horizon to infinite 
into (— 00 , 00 ) region, and r* is defined as 


dr* = / 3dr . (74) 

The effective potential [29] l3U]is 

V(r) = ^[1(1 + 1) + c/( 1 - s 2 )r + (s - s 2 )( 1 - a)], (75) 

where s = 0 and 1 corresponds to the scalar field and Maxwell field, respectively. Here, l 
is the multipole quantum number (the angular harmonic index). The effective potential 
has a maximum at r* = (r*) 0 . 

The WKB approximative method was first applied by Schutz and Will [31] to the 
problem of scattering around black holes. This method is based on matching of the 
asymptotic WKB solutions at spatial infinity and the event horizon with the Taylor 
expansion near the top of the effective potential barrier through the two turning points. 
Here, we adopt the third-order WKB method [28]. The QNMs is as follows 

w 2 = [Vb + (—2K 0 ")U] - iu{- 2V r 0 ")3(l + ft), (76) 


where the subscript 0 on a variable denotes the value of the variable at (r *) 0 {Vff 7 ^ 0 ), 
and v = n + -7 (overtone number n = 0,1, 2, • • ■, for Re{uj ) > 0). Furthermore, 



(_ 2 t / 0 ") 1 / 2 16912 VRq 



(19 + 28 z/ 2 ) 


( 78 ) 
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Table 1. QNMs of Type I A furry black hole for electromagnetic perturbation with 


r s 

> 0, S = Ar); and A 

, = 2. Here, we 

have worked in unite with |r s | = 1. 

l 

n = 0 

n = 1 

n = 2 

n — 3 

1 

0.3782-0.1643* 

0.3007-0.5316* 

0.1935-0.9198* 

0.0480-1.3174i 

2 

0.7207-0.1685* 

0.6725-0.5190* 

0.5970-0.8888* 

0.5010-1.2685* 

3 

1.0404-0.1698* 

0.0059-0.5161* 

0.9466-0.8758* 

0.8704-1.2463* 

4 

1.3537-0.1703* 

1.3269-0.5150* 

1.2784-0.8690* 

1.2141-1.2327* 

5 

1.6644-0.1706* 

1.6424-0.5145* 

1.6015-0.8651* 

1.5458-1.2238i 


Table 2. QNMs of Type I A furry black hole for electromagnetic perturbation with 
r s >0, S = Ts and A = Here, we have worked in unite with |r s | = 1. 


l 

n = 0 

n = 1 

n = 2 

n = 3 

1 

0.1656-0.0507* 

0.1512-0.1579* 

0.1310-0.2703* 

0.1050-0.3844* 

2 

0.3000-0.0514* 

0.2913-0.1563* 

0.2769-0.2648* 

0.2590-0.3759* 

3 

0.4287-0.0516* 

0.4225-0.1559* 

0.4113-0.2624* 

0.3968-0.3714* 

4 

0.5557-0.0517* 

0.5508-0.1557* 

0.5418-0.2613* 

0.5296-0.3687* 

5 

0.6819-0.0517* 

0.6780-0.1556* 

0.6704-0.2606* 

0.6599-0.3671* 


where the primes and the superscript (j) denote differentiation with respect to r*. The 
numerical results of QNMs for electromagnetic perturbation (s = 1) are listed in Tables 
1-6. As a reminder, the results of QNMs verifies our conjecture: Type I is stable and 
Type II is unstable. 

It is worthy that the corollary of QNMs are the same as those for the massless s = 2 
perturbation. However, we will run up against difficulties if we use WKB method to 
the scalar perturbation (s = 0). In Type II and s = 0 case, the effective potential (175]) 
is not a potential barrier for lower l so that it can’t reduce to the problem of scattering 
around black hole. 

For Type I furry black holes, the real parts of the quasinormal frequencies increase 
as l, which means that the large l is, the faster black hole oscillates; the imaginary parts 
are always negative, which corresponds to the oscillation of the black hole decays. Thus, 
Type I black hole will tend to stable eventually. For Type II furry black holes, there exist 
the positive imaginary parts for the lower modes so that they are unstable. In Tables 1-8, 
we have taken \r s \ = 1, therefore, the unite of QNMs is |r s | _1 ~ 1.02 x 10 5 (M Q /M)sec _1 . 
For ( l,n ) = (1,0) quasinormal mode of Type IIA furry black hole with \M\ = M 0 
in Table 5, the electromagnetic perturbation will enhance 42 times per microsecond. 
Obviously, the life of unstable black hole is proportional to its mass. 

5. Phenomenological consequences 

5.1. A new astronomical object and its behavior 

The behavior of the metric (]63l) is determined by the two integral constants r s and Uo , 
and the values of parameters Co and C\. At the origin r = 0 the term proportional to r s 
is singular, so the metric always possesses a singularity unless r s = 0. If we take r s = 0, 
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Table 3. QNMs of Type IB furry black hole for electromagnetic perturbation with 


r s 

> 0, S = — and A = 2. Here, we 

have worked in 

unite with |r s | = 1 

l 

n = 0 

n = 1 

n = 2 

n = 3 

1 

0.6082-0.1867? 

0.5541-0.5836? 

0.4796-1.0012? 

0.3851-1.4244? 

2 

1.1064-0.1908? 

1.0745-0.5798? 

1.0222-0.9825? 

0.9571-1.3945? 

3 

1.5819-0.1918? 

1.5594-0.5792? 

1.5192-0.9750? 

1.4669-1.3793? 

4 

2.0510-0.1923? 

2.0336-0.5790? 

2.0012-0.9713? 

1.9573-1.3704? 

5 

2.5172-0.1925? 

2.5030-0.5789? 

2.4760-0.9693? 

2.4384-1.3648? 


Table 4. QNMs of Type IC furry black hole for electromagnetic perturbation 
r s < 0, S = §(— r s )z and A = Here, we have worked in unite with |r s | = 1. 

1 

n = 0 

n— 1 

n = 2 

n = 3 

1 

0.0876-0.0196? 

0.0843-0.0599? 

0.0791-0.1017? 

0.0938-0.1440? 

2 

0.1552-0.0197? 

0.1533-0.0596? 

0.1499-0.1003? 

0.1457-0.1418? 

3 

0.2207-0.0198? 

0.2193-0.0595? 

0.2167-0.0998? 

0.2133-0.1406? 

4 

0.2207-0.0198? 

0.2844-0.0595? 

0.2824-0.0995? 

0.2795-0.1400? 

5 

0.3501-0.0198? 

0.3492-0.0595? 

0.3475-0.0994? 

0.3451-0.1396? 


Table 5. QNMs of Type IIA furry black hole for electromagnetic perturbation with 
r s < 0, S = and A = 2. Here, we have worked in unite with |r s | = 1. 


1 

n = 0 

n = 1 

n = 2 

n = 3 

1 

13.243+36.813? 

41.621+98.291? 

90.522+172.11? 

156.37+257.15? 

2 

9.4002-24.533? 

9.6272+85.761? 

42.807+150.31? 

88.781+221.11? 

3 

31.462-22.446? 

13.037-80.206? 

12.991+141.17? 

48.945+205.92? 

4 

48.402-22.892? 

32.260-77.273? 

10.168-136.13? 

19.633+197.64? 

5 

63.455-23.280? 

49.602-75.662? 

30.041-132.91? 

4.3625-192.43? 


Table 6. QNMs of Type IIB furry black hole for electromagnetic perturbation with 
r s > 0, S = — lOrf and A = Here, we have worked in unite with |r s | = 1. 


1 

n = 0 

n = 1 

n = 2 

n = 3 

1 

1199.1+1913.6? 

2779.6+4486.8? 

5128.1+7685.5? 

8234.5+11451? 

2 

591.93+1512.2? 

1551.1+3811.0? 

3161.7+6517.5? 

5367.0+9582.8? 

3 

143.45+1235.1? 

801.08+3511.2? 

2042.4+6018.4? 

3783.2+8769.7? 

4 

327.47-1018.3? 

229.45+3341.3? 

1247.2+5747.2? 

2691.9+8322.6? 

5 

809.10-930.26? 

252.12-3233.6? 

618.35+5580. 2i 

1855.5+8046.3? 


Table 7. QNMs of MBH with the metric (RI51) . e = yA. anc j Here, we have 


worked in unite with |r s | = 1. 


1 

n = 0 

n = 1 

n = 2 

n = 3 

1 

0.4917-0.1862? 

0.4226-0.5916? 

0.3286-1.0180? 

0.2038-1.4509? 

2 

0.9141-0.1901? 

0.8715-0.5818? 

0.8045-0.9915? 

0.7209-1.4110? 

3 

1.3133-0.1912? 

1.2828-0.5795? 

1.2300-0.9799? 

1.1627-1.3908i 

4 

1.7058-0.1917? 

1.6820-0.5786? 

1.6389-0.9738? 

1.5817-1.3782? 

5 

2.0954-0.1919? 

2.0760-0.5781? 

2.0396-0.9703? 

1.9900-1.3699? 
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Table 8. QNMs of MBH with the metric (l63l) . e = and ^ Here, we have 

worked in unite with |r s | = 1. 


l 

n = 0 

n= 1 

n = 2 

n = 3 

1 

0.4917-0.1862* 

0.4226-0.5916* 

0.3286-1.0180* 

0.2039-1.4509* 

2 

0.9142-0.1901* 

0.8716-0.5818* 

0.8046-0.9915* 

0.7210-1.4110* 

3 

1.3134-0.1912* 

1.2828-0.5795* 

1.2301-0.9799* 

1.1628-1.3908* 

4 

1.7059-0.1917* 

1.6821-0.5786* 

1.6390-0.9738* 

1.5818-1.3782* 

5 

2.0956-0.1919* 

2.0761-0.5781* 

2.0398-0.9703* 

1.9901-1.3700* 


1.5 


1.0 


0.5 


0.0 

0 1 2 3 4 5 

Log r 


Figure 1(a) 



the metric tends to 1 since the logarithmic term becomes m 2 /2vf at r = 0. In the case 
of r 1, a tends to 1 — {m 2 /2vff) which describes as a metric with a solid angular 
deficit [32] [33, |33J [35] since the logarithmic term is r~ 2 order as r —* oo. 

There is another singularity at r = a/uo for the metric field a. The point r = a/uo 
is a pole of order 2 for the scalar curvature R. These singularities may or may not be 
hidden by the horizon depending on the values of r s and Uq. The solutions possessing 
the horizon are candidates for modified black holes (MBH). 

The horizon is always present if the integral constant r s > r crit for fixed Uq, a and 
m (see Figure 1), where 


a V2e 3 u 0 

Tcrit 1 d (■h 

Uq (1 — £ Z )m 


(79) 


Actually, it is easy to prove the above-mentioned result by the asymptotical behavior 
of gravitational held. 

Next, let us discuss stability of MBH with the metric (j63|) . Using (163|) and (1751) - 
), we list the results of QNMs in the Tables 7 and 8, which show this MBH is stable. 


5.2. The deflection of a light ray and the perihelion shift of Mercury 

For a given gravitational held of the astronomical object with the metric (1631) . the 
equation of motion then permits us to predict the trajectories of light signals, which we 
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Figure 1(b) 



Figure 1(c) 



Log r 


Figure 1. The existence of horizon is examined for the metric (j63j) . Figure (a) 
corresponds to naked singularity at r = a/uo with the integral constant r s < r crit . The 
horizon rh is shown in the figures (b) and (c) where the integral constant r s = r C rit 
(b) and r s > r cr i t (c). 


can compare with observation data. In the case of e -C 1, we should be able to neglect 
the e 3 term in (j63l) then the geodesic equations in the plane 6 = n/2 are 

(80) 

(81) 

(82) 


dt 1 — £ 2 
dp a(r ) ’ 

2 # = 7 

dp 

^(< )2 


1 - £ 2 J 2 - 
12 _ ^ , 


a(r) m 

where p is an affine parameter to coincide with the coordinate t asymptotically, and 
the integral constant E is such that ds 2 = —E 2 dp 2 . J is also a constant of motion, the 
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angular momentum. By rescaling the following quantities, 

r s = r s (l - £ 2 ), 0 = 0(1 - £ 2 )~K J = J( 1 - £ 2 )K 

E 2 = E' 2 ( 1 - e 2 ), t = t( 1 - £ 2 )-\p = p( 1 - e 2 )~\ (83) 

m-m l°°k same as those in the ordinary Schwarzschild metric. Therefore, we 
can use results of the Schwarzschild case for our issues. In a Schwarzschild metric, the 
change in the angular coordinate for a light ray scattered is A0 = 7r + 2r s r0 1 , where ro 

is the closest distance of the light ray from the centre of gravitational field. The angle 

by which light is deflected for this gravitational field is 

A0 = tt (1 - e 2 y 1/2 + 2r s (l - e 2 )" 372 ^ 1 - i r. (84) 


One of the classical Einstein effects is the deflection of a light ray in the gravitational 
held of the sun. The effect is small but large enough to the detectable, which 
has been tested by observations during total eclipses of the sun on the apparent 
positions of stars whose light has passed close to the limb of the sum. The result 
of optical measurements is consistent with the prediction of GR, but hardly a precision 
conformation. However, the optical measurement can be substituted with the radio 
interferometry in these measurements, which has led to the development of very long 
baseline interferometry (VLBI). Using radio transmission from certain quasars and 
measuring the deflection as the source is eclipsed by the sun, the most precise result is 
that A(f exp /Af>GR = 1-0001 ± 0.0001 in terms of a ratio of the observed deflection and 
the theoretical value A cfoR = 4 GM Q /R Q ~ 8.488 x 10 -6 rad for GR [36] . Consider the 
case of r 0 = R Q in (l84]h we can write e 2 in terms A <f G R and A <f> exp as follows 


= !-{[y 


7T 


A0 


exp 


A(pGR A<Pgr 


)+VI]i + [- 


7T 


A0, 


exp > 


2 A <f>GR A<f>GR 


VA]^}- 2 ,(85) 


where 


A = 


7T 


A 0, 


exp \ 2 


r + ^(- 


7T 


( 86 ) 


4 V A (pGR A(pGR 27 A (f>GR 

Using (T53T) and Af> exp < l.OOO2A0 GR , we have the upper bound of the parameter e 2 , 
that is, e 2 < 1.081 x 10 -9 . 

There are other constraints from the solar system tests such as the perihelion shift 
of Mercury and the gravitational redshift. For the astronomical object with the metric 
), the perihelion shift 


A0 = 


67 tGM 


(87) 


ac 2 (l — e 2 )(l -e 2 ) 2 ' 

where e is the eccentricity of the ellipse and a is semi-major axis. On the other hand, the 
measured perihelion shift of Mercury has been known accurately. After the perturbing 
effects of the other planets have been accounted for, the excess Af>GR is known to about 
1CU 3 from radar observations of Mercury between 1966 and 1990 [3?j|. Therefore, we 
have the upper bound of the parameter e 2 , that is, £ 2 < 1.498 x 10^ 3 . This constrain is 
weaker than one from the deflection of a light ray. 
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The geometric interpretation of e 2 is a solid deficit angle of space. It is useful 
that the constraint is converted into a constraint on a more physically transparent 
parameter. We note that the dimension of uo is [Length]. In the M p i and A terms, we 
have Uq > 25.84A 2 M~ l 1 from the solar system tests. The existence of black hole is the 
integral constant r s > r cr i t , which can be rewritten as 

eM P i 

-o 


_ / Ci l 

“ 2cV 4 (1 -e 2 )A 2 ’ 


( 88 ) 


by using the parameters in the Lagrangian (JT]) such as A, cq and C\. 


5.3. Einstein ring 


If the source of light, the astronomical object with the metric fib.311 . and the observer 
are aligned exactly, all the rays that pass at the appropriate impact parameter around 
this object, at any azimuth, reach the position of the observer. Under these special 
circumstance, the observer sees an infinite number of images, which form a ring around 
this object. Assuming the source is much farther from this object than the observer, 
the rays incident on this object are then nearly parallel to the line of the alignment, and 
the deflection angle required for the ray to reach the observer is b/D, where D is the 
distance from this object to the observer. Thus, the angular radius of the Einstein ring 
is 

b £ 2 tt 2 r, 

D = ~2~ + IT’ 

which can also be written as 

b £ 2 7 t je A i r 2 2r~ s 

D~ "IfT + 1 T' 

If e 2 = O(10 -9 ), we obtain 




b 

15 


> O.larcsec. 


(91) 


In the Schwarzschild case, the deflection angle tends to zero as D tends to infinity. In 
our case, the deflection angle tends to £ 2 7 t/2 as D tends to infinity, that is to say, b also 
tends to infinity with D because the solid deficit angle extends to infinity. This angle 
is at the limit of resolution of optical telescopes, but it is well within the resolution 
attainable by radio telescopes. Such an Einstein ring should be observable with ratio 
telescopes if they exist. 


5-4- Cosmological density constraint 

In order to obtain a cosmological constraint on the abundance of this astronomical 
object, we use the upper bound of the density parameter 0 = p s /p cr it, o- Larson and 
Hiscock have estimated the astrophysics bounds on cosmic strings by the contribution of 
energy density [38] . Similarly, we can estimate the minimum separation between these 
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astronomical objects. If the typical distance between these objects is defined to 2L, the 
average density np to the 0(e 2 ) order is 

1 r s + e 2 (L - R s ), 


ps ~2G [ 




(92) 


4vrL 3 /3 

where R s is the radius of this astronomical object. From p Sj0 < Per up and r s , R s -C L, 
we have 


L > snp 0 Hp. 


(93) 


where Hq is Hubble parameter at present. Obviously, the separation for these objects 
is more larger than the size of the Galaxy, so it is not easy to find these objects if they 
exist unless e 2 is tiny (e 2 < 10~ 12 ). 


6. Conclusion and discussion 

In this paper we have developed a detailed study of the spherically symmetric solutions 
in Lorentz breaking massive gravity. We have shown clearly that all solutions of 
cLF = 0 form a commutative ring R T if the Stiickelberg held <f l is taken as a hedgehog 
configuration 0* = f>(r)x l /r. Under the hedgehog configuration and J- is a polynomial, 
there is the only solution <f = br for Einstein equations. Moreover, we have obtained the 
expression of solution to the functional differential equation with spherically symmetry 
if T G R T . Using this universal formula, we give known solutions including the 
Schwarzschild, dS and AdS if we take T G S T where S T C R T and dTjdX = 0. 
When T e R T but T 0 S T , we have given some analytical examples including the furry 
black holes and a new metric solution fl63jl . If we take r s > r cri t in the metric (l63]h this 
solution will describe a modified black hole or star. The stability of the Stiickelberg 
held 0* is guaranteed by the topological one. The stability of these black holes under 
perturbations are discussed using the analysis of Komar integral and the results of 
QNMs. We also discussed some phenomenological consequences for these solutions. 

We have discussed the minimum separation between the astronomical objects from 
the cosmological restriction in detail. It shows that if the typical distance between these 
objects is defined to be 2 L, L > Obviously, it is not easy to find these objects 

because the separation for these objects is more larger than the size of the Galaxy. 

Furthermore, one can consider a region where matter is present in the form of 
a perfect fluid with a constant energy density p and pressure p. In other words, the 
novel features can be ascribed to the Stiickelberg fluid turned on by matter inside the 
body. Thus, a self-gravitating body can be described by matching the exterior with 
the interior solution. In comparison to GR, the resulting equations are difficult to solve 
analytically even for a constant matter fluid because of the presence of the Stiickelberg 
fields. Therefore, one has to rely upon a nonstandard perturbation expansion ra- It is 
an interesting question whether there is a star solution with pure Stiickelberg held. We 
will study this question further and it is hopeful that many observations will be done 
to test LB massive gravity. 
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Appendix A. Stiickleberg fields and subring 


The functionals JC(T) and 'H(J-) are defined as 

/C(jr) = A - 2AJ7, (A.l) 

H(T) = + (A. 2) 

ry* A ry->et ry*\J 

then the Einstein equations can be rewritten as 

a' a — 1 m 2 . _ , . . 

~ + ~2 - yK-iF) = 0, (A.3) 

III 9 

a ot m* 

Y + 7 - —nr) = o, (A.4) 

0' 2 20' 4 30' 6 

XTx + yJ 7 i- yJ 7 2 + yJ 7 3 — 0. (A.5) 

In this appendix, we shall prove the following lemmata using (lA.3MA.5jh which describe 
substructure of Vf and the Stiickleberg fields. 

Lemma 1. Let T = Y^ a hi 1 i 2 h{^) l ° wl i wl 2 wl i i s a polynomial (Z 0 , , Z 2 , Z 3 > 0) and 
T G Vf ", then <f l = bx l is the only solution for the spatial Stiickleberg fields under 
the static spherically symmetric ansatz ([3]). The solution T can be classified by the 
parameter s = l\ + 2 l 2 + 3Z 3 . 

Proof The solution T is generally a function of X, Wi and 0 /2 for flA.5|h that is 
to say, T = A(A, wi, w 2 , w 3 , 0' 2 ). However, T = T{fX,w\,w 2 ,w 3 ) is necessary for the 
SO( 3) symmetry in the 0* internal space [11] . Thus, 0* = bx l is the only solution for the 
spatial Stiickleberg fields under ansatz (J3]). In this case, (IA.5D is reduced to an algebraic 
equation as follows 

^(-l) Zl+ ' 3 a Wl y 3 fe 2 s (^)' 0 +s M(A)(l+2A) /l (l+2A 2 )' 2 (l+2A 3 ) Z3 = 0(A.6) 


where 


M{ X) = l 0 + 


'A 

1 1 2/2 


+ 


+ 


3/.? 


(A.7) 


1 + 2A 1 + 2 A 2 1 + 2 A 3 ’ 
and s = l\ + 2/ 2 + 3 / 3 . Therefore, the polynomial solution T can be classified by 6 2s .D 
It is easy to find the linear independent T for s = 0,1, 2, 3 as follows 

s = 0 : A = C 0 1 ; (A.8) 

b 2 

s = 1 : T = Cl(— + Wi ) + b 2 x (s = 0 solution); (A.9) 

JC 

s = 2 : T = C\{wl — W 2 + 4:b 2 wi) + b 2 x (s = 1 solution); (A. 10) 
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s = 3 : 


7 

7 

7 

7 

7 


^2(^0 — w z) + b 2 x (s = 1 solution); (A. 11) 

A 2 

C\(w\ — 3 w\W 2 — 6 b 4 wi + 2 W 3 ) + b 2 x (s = 2 solution); 


(A. 12) 

b e 

Cft (+ w 3 ) + b 2 x (s = 2 solution); (A. 13) 

b 6 2b 6 2 6 6 

+ ^1^2) + b 2 x (s = 2 solution); 

(A.14) 

C 3 (— + — + — + wf) + b x (s = 2 solution). 

(A.15) 


where Cj are constants. Using addition and multiplication of ring, we can obtain solution 
7 for s > 4. 

The argument is given in the Lemma 1 that 7 = 7(X, wi,W 2 ,Ws) is necessary for 
50(3) symmetry in the 0 l internal space. In other words, 7 cannot explicitly depend on 
<f>'. Although the quantities themselves depend on ft (see (USD), 1(11,2,3 are reduced 
to the functions of A -1 iff 0 = hr. 

A reasonable theory of gravity ought to contain the known spherically symmetric 
solutions in GR (Schwarzschild, Sch-AdS or Sch-dS solution). From (1A.3I) and (1A.4I) . 
we have 


a 

~2 


«-!_ 2 (YT ^ T 2 < i ^ 4 77 , 77 \ 

— rn (X7x + —7 1 — —7 2 + 3 )- 


r 


(A.16) 


iy> Zj ^ 

Using (1A.5D . (IA. 16[) can be written in the form 

“- 1 - 9 > 4 <Ps T WA17^ 

—3— - m [(— - —)7 1 - 2(— - —^)7 2 + 3(— - -t7z)7 3].(A.17) 


a 

~2 


X 


A 2 ' 


A 3 ' 


In the case of 0 = br and A = 1, (1A.17|) is reduced to 


and 



(A.18) 


a = 1 — — + Ar 2 , (A.19) 

r 

where r s and A are integral constants. 

In the aforesaid considerations, we have used the requirement of vacuum solutions 
in GR. Conversely, we should consider what is restriction of 7 if there is only 
the vacuum solution. We find that the constraint is simple 7(7) = 7(7), and 
S T = {7 £ Vf | 7(7) = 7(7)} forms a subring, S T C Vf. Here, we have chosen 

F = E a ioiihi3(jt) l0w i w 2 w 3 ( l o, h, h, h > 0). 

Lemma 2. If we take 7 £ S the Einstein equations (IA.3D - (IA.5[) have only 
the Schwarzschild, Sch-AdS or Sch-dS solution. S T has an equivalence definition 
S T = {7 £ R F \7x = 0}. 
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Proof. Combining (I A. 31) and (1A.4j> with IC(F) = T-L(F), we find that a satisfies 
(1A.18I) . Therefore, the Einstein equations have only the Schwarzschild, Sch-AdS or Sch- 
dS solution for any F G S T . Thus, X must be 1 and l-L(F) = iC(F) is equivalent to 
Tx = 0 from the Einstein equations. □ 

Lemma 3. Let F = Y^ a ioiR*h{jc) l ° wl i wl 2 wl i is a polynomial (lo,h,l 2 ,h > 0) 
and S T = {F G Vf \ JC(F) = f-i(F)}, then S T forms a subring of Vf. The set 
Sff = {F G TZ T | IC(F) = cH{F),c A 1 is a fixed constant} forms a subgroup of the 
abelian group . 

Proof. Suppose that /\/ 2 G S T and K (/*) = %(/*), (i = 1,2). From Lemma 2 
and (IA.5D . we have /Cf]/ 1 + p) = Pip 1 + / 2 ) and /C(/ 1 / 2 ) = 'H(/ 1 / 2 ) so S' T C Vf. For 
f 1 , p G S^f, we only have /C(_/‘ 1 + f 2 ) = cH(P + / 2 ) so S)f C 

Note that we always obtain new solutions beyond the furry black holes using the 
addition and multiplication for F 1 G Sfi and F 2 G SF 
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